Order by distortion and string modes in pyrochlore antiferromagnets 
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We study the effects of magnetoelastic couplings on pyrochlore antiferromagnets. We employ 
Landau theory, extending an investigation begun by Yamashita and Ueda for the case of S = 1, and 
semiclassical analyses to argue that such couplings generate bond order via a spin-Peierls transition. 
This is followed by, or concurrent with, a transition into one of several possible low-temperature 
Neel phases, with most simply collinear, but also coplanar or mixed spin patterns. In a coUinear 
Neel phase, a dispersionless string-like magnon mode dominates the resulting excitation spectrum, 
providing a distinctive signature of the parent geometrically frustrated state. We comment on the 
experimental situation. 



Geometrically frustrated magnets (l]-||] are examples of 
strongly interacting systems: the vast degeneracy of their 
classical ground states makes them highly susceptible 
even to small perturbations. By analogy with quantum 
Hall systems, where the Landau levels are also macro- 
scopically degenerate, one expects a variety of phases in 
perturbed frustrated magnets, from Neel states to spin 
glasses or liquids, with valence-bond solids along the way. 

Probably the world's most frustrated spin system is the 
classical Heisenberg antifcrromagnet on the pyrochlore 
lattice (Fig. |l|) where spins reside at vertices of tetrahe- 
dra. The number of its classical ground states, which 
are attained when total spin on each tetrahedron Stot = 
^^^-^ Si = 0, is so large that, exceptionally, it does 
not order at any finite temperature Q|. In real com- 
pounds, deviations from the classical Heisenberg model 
(e. g. dipolar interactions, single-ion anisotropy or quan- 
tum fluctuations) determine which ground state is se- 
lected at the lowest temperatures. 




FIG. 1. The pyrochlore lattice 

In this note, we discuss an elegant mechanism for lift- 
ing the frustration through a coupling between spin and 
lattice degrees of freedom. The high symmetry of the 
pyrochlore lattice and the spin degeneracy drive a dis- 
tortion of tetrahedra via a magnetic Jahn-Teller ( "spin- 
Teller") effect. The resulting state exhibits a reduction 
from cubic to tetragonal symmetry and the development 
of bond order in the spin system with unequal spin cor- 



relations (Si • Sj) on different bonds of a tetrahedron. In 
the ordered phase, there are 4 strong and 2 weak bonds 
per tetrahedron — or vice versa. This phenomenon was 
uncovered by Yamashita and Ueda for pyrochlore an- 
tiferromagnets with spins 5 = 1, for which they described 
an AKLT-style wavefunction with the requisite bond 
order. In the following we study this phenomenon in 
the semiclassical limit with added insight from Landau 
theory, and discuss its consequences for the excitation 
spectrum. As many of the candidate systems have mod- 
erately large spins and order at finite temperatures, our 
methods should work well — in particular, they allow us to 
treat the Neel order that can (and experimentally does) 
appear in addition to the bond order. 

We begin by identifying a two-component bond order 
parameter at the level of a single tetrahedron and pro- 
ceed to construct its Landau free energy for states of the 
infinite lattice. In one family of cases the spin-Peierls 
transition can be first order, as indeed observed in many 
pyrochlore antiferromagnets |7|-p^. 

Once the system is in a bond-ordered phase, the low- 
ered symmetry of valence bonds reduces frustration of 
spins and sets the stage for Neel (spin) ordering. In gen- 
eral, we expect a separate phase transition into an anti- 
ferromagnetic state. However, the first-order spin-Peierls 
transition may turn the system directly into an antifcrro- 
magnet, bypassing a metastable spin-Peierls phase. Both 
possibilities seem to be realized in different materials 
(|ll|,|2| . The particular type of the Neel order depends on 
the underlying bond ordered state. While collinear order 
is most easily generated, we show that unusual coplanar 
states are also possible. 

Finally, we obtain a striking signature of the Jahn- 
Tcller distortion in collinear Neel states. The spectrum 
of spin excitations in the distorted antifcrromagnet con- 
tains a large number of modes clustered near a finite 
frequency. These magnons, a remnant of pyrochlore zero 
modes, are confined to strings of parallel spins. We there- 
fore call them string modes. A resonance highly reminis- 
cent of such modes has been observed by S.-H. Lee et al. 
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p2[ in ZnCr2 04. The strings should also be observable 
in ZnV2 04 where a structural distortion and a coUinear 
Neel order at low temperatures have been firmly estab- 
lished . A weaker feature might exist in the itinerant 
magnet YMn2. We turn now to the details of our results. 

Single tetrahedron. Yamashita and Ueda |^ have 
discussed the Jahn-Teller distortion on a single tetrahe- 
dron for the case of spins S = 1/2. Their analysis can 
be applied to general spin virtually unchanged. On an 
undistorted tetrahedron, all exchange constants are equal 
and any state with total spin Stot = minimizes the clas- 
sical Heisenberg energy J X)i> j ' ■ To see which of 
these is selected in the presence of a small spontaneous 
distortion, we minimise the sum of spin and elastic en- 
ergies. The energy of a bond JijSi ■ Sj depends on ionic 
displacements through variation of Heisenberg exchange 
Jij with the length of interatomic bonds and bond an- 
gles. Expanding magnetic and elastic energies to lowest 
order in the displacements {a = 1 . . . 12) we obtain 

E = ^{dJij/dXa){Si-Sj) Xa+^ ka(3XaXp/2. (1) 

As the reduction of magnetic energy is linear in Xa , it will 
beat the quadratic cost in elastic energy. Hence a Jahn- 
Teller distortion will result from the spin degeneracy. 




FIG. 2. Possible values of the two components of the 
bond vector f = (/i,/2) = {fcosO, fsinO) are bounded by 
an equilateral triangle in the (/i, /2) plane. Also shown are six 
extremal spin configurations (filled and open circles). Strong 
(weak) bonds are denoted by solid (dashed) lines. 

For tetrahedral (T^) symmetry, displacements Xa that 
contribute to the magnetoelastic energy (^) must belong 
to the same irreducible representation as the six bond 
variables Si ■ Sj . These bond variables can be thought of 
as classical quantities or, for quantum systems, as expec- 
tation values of the corresponding operators; our symme- 
try analysis applies to both cases. 



In the experimentally relevant regime well below the 
Weiss temperature, Stot ~ 0, and we are left with only 
three independent quantities. One of them, the ground- 
state energy X]i> j ' ^ invariant under all symme- 
try operations, induces a uniform rescaling of the tetra- 
hedron (which is not a distortion). The other two, 
(/ii /2) = f, measure the disparity of the bond variables: 

/i = [(Si + S2) • (S3 + S4) - 2Si • S2 - 2S3 • S4] /VT2, 
/2 = (Si • S3 + S2 • S4 - S2 • S3 - Si • S4)/2. (2) 

They transform as the E representation of the group Td 
and couple to phonons of the same symmetry. For a given 
configuration of spins, the minimum magnetoelastic en- 
ergy (m) can be written in two ways (C^ = const): 



E ~idJ/dxE)^f/{2kE) ^C^^Ci ^(S, 



(3) 



here dJ/dxE and are the appropriate magnetic 
and elastic constants. Minimization of total energy is 
achieved when the vector f has maximum length. From a 
semiclassical analysis, we find that f is restricted to lie in 
an equilateral triangle, at the corners of which its length 
is maximized (filled circles in Fig. |^). This corresponds 
to two weakened and four strengthened bonds, a result 
which is easily rationalised in terms of classical ground 
states of the tetrahedron: only one of them (modulo sym- 
metries) fully satisfies four bonds and completely frus- 
trates the remaining two, which are then, respectively, 
strengthened and weakened by the distortion. This state 
is the coUinear state, which is of course favoured by the 
effective "biquadratic exchange" (||). 

Infinite pyrochlore lattice. The Jahn-Teller distortion 
of individual tetrahedra will analogously drive a spin- 
Peierls phase transition on the infinite pyrochlore lattice. 
There now exist an infinite number of phonons that could 
become soft. We here discuss the simplest but already 
very rich case of a phonon condensate with lattice mo- 
mentum q = 0, for which all tetrahedra of a given type 
distort in the same way. The two inequivalent tetrahe- 
dra (type A and B) which make up the pyrochlore lattice 
reside on a bipartite (diamond) lattice (Fig. |l|). 

The symmetry group of the pyrochlore lattice, Oh, is 
enlarged from by the operation of inversion through a 
site, which exchanges A and B tetrahedra: Td®Ci = Oh- 
Irreducible representations of the cubic group Oh are 
those of Td with an additional quantum number, par- 
ity under the inversion. Thus, there are two doubly de- 
generate displacement modes that couple to bonds in a 
nontrivial way. The even phonon doublet Eg creates a 
uniform distortion of the entire lattice; the odd doublet 
E^ causes a staggered distortion (f^ — — fs). 

Generalising Eq. ^ for the energy, using small sub- 
scripts for parity and capitals for tetrahedron type, gives 
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For kg = ku, fA and /b are separately maximized. For a 
stiffer odd phonon, ku > kg, fA and fe are in the same 
corner of the triangle in Fig. ||. When kg slightly exceeds 
ku, "antiferromagnetic" coupling puts and fs in two 
different corners. Fig. |^(a,b) displays spin arrangements 
corresponding to these two types of bond order. The 
latter has been observed in YMna and MgV204 |0,|,|lO|- 
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FIG. 3. (a-b) Sample coUinear states of classical spins 
on the pyrochlore lattice, projected along (001). Frustrated 
bonds are shown as broken lines. They form (a) infinite chains 
and (b) infinite spirals, (c) Neel order with coUinear spins on 
tetrahedra A and coplanar spins on tetrahedra B. 

We next discuss nature of spin-Peierls phases and tran- 
sitions by considering the Landau free energy A{g,u), 
where g, u are the uniform and staggered bond order pa- 
rameters jl^ , g, u = (fyi ± £5)72. In a paramagnetic 
state, the bond order parameters g, u = 0; note that 
spin-Peierls order g, u ^ does not necessarily imply a 
Neel order. The terms allowed by Oh symmetry are 

^(g, u) = Ggg'^ + bgg^ cos 39g + Cgg^ + ... 

+ auU^ + CuU^ + duU^ cos 66'„ + e„u^ + . . . (5) 
+ buU^gcos{2eu + 9g) + . . . 

Quite generally, one can read off that the phase transition 
will be discontinuous if it is driven by the even phonon 
Eg: the symmetric cube [E'^] contains the trivial repre- 
sentation making it possible to write a cubic invariant 
~ 352) = g^ cosWg in polar coordinates of Fig. ||. 
On the other hand, the transition will be continuous if 
it is driven mainly by the odd phonon We also note 
that both g and u can be present owing to a nonlinear 
coupHng. Then 7^ f_A 7^ fs 7^ 0, as in the case depicted 
in Fig. |^(b). Below the transition temperature, either Og 
or a„ becomes negative. Omission of higher-order terms 
requires that Cg, c„ > 0. Signs of bg,bu, and d^ will deter- 
mine the nature of spin-Peierls phases. In our classical 
model, bu > 0, which we assume in what follows. 

For the phase transition driven by the uniform distor- 
tion Eg, only the first line of Eq. ^ needs to be taken 
into account. Depending on the sign of bg, the ordered 
phase features either two strong and four weak bonds per 
tetrahedron [Fig. |^(a)], or two weak and four strong ones 
[Fig. ^(b)]. In the ordered phase, a growing g softens 
the staggered mode u through mode coupling, the third 



line in Eq. |[ Once g > a„/6„, there is a second transi- 
tion into a phase with both g, u 7^ (f^ 7^ fs)- The two 
possible sequences of transitions are shown in Fig. ^(c,d). 

When the staggered phonon Eu drives the transition, 
g can be neglected initially, we are then left with the sec- 
ond line in Eq. |5[ This transition is continuous (unless 
the g mode is nearly unstable) . The sign of d„ determines 
the direction of the order parameter u. The distinct pos- 
sibilities are 9u — (two strong bonds on tetrahedra A 
and four strong bonds on tetrahedra B) and 9u = 7r/2. 
The free energy in the ordered state with u 7^ shows 
that a subdominant order parameter g of 0{u^) arises 
simultaneously with u [Fig. ^(e,f)]. 




FIG. 4. Evolution of the spin-Peierls order parameters 
fA ~ g -I- u and fs — g — u. (a-b) g condenses, u = 0. 
(c-d) g condenses first, u appears later, (e-f) u condenses, 
g is generated as a subdominant order parameter. In the 
paramagnetic phase (filled circle), = fs = 0. 



Neel order. Lowering the bond symmetry relieves frus- 
tration and thus encourages spin ordering [Q. However, 
as bond and spin ordering are not required to coincide 
even classically, the former will likely precede the latter. 
Indeed, frustrated systems are the ideal setting for this to 
happen: the classical ground state degeneracy supplies a 
plethora of local zero-energy modes [|| , which can allow 
spins to fluctuate strongly even in the presence of bond 
order. This, essentialy classical, mechanism is distinct 
from the more conventional avenue of obtaining valence- 
bond solids through enhanced quantum fluctuations for 
small spins. Nonetheless, if the spin-Peierls transition is 
discontinuous, it can take the paramagnet directly into 
the Neel phase, past the metastable spin-Peierls state. 
Both variants seem to occur in nature ||ff|,Q . 

The phase transition from a bond-ordered state to a 
Neel phase can also be analysed within the Landau frame- 
work by adding spin averages (Si). For the sake of sim- 
plicity, we make the assumption that the contribution 
of the spin condensate (Si) • (Sj) will tend to reinforce, 
rather than weaken or alter, bond order (Si • Sj). Then 
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one can use Fig. || as a guide to identify the appropri- 
ate low-temperature Neel phase for each of the spin- 
Peierls phases of Fig. ^. The spin-Peierls states shown 
in Fig. l4|(a,d) exhibit the coUinear Neel order shown 
in Fig. 0(a,b), respectively. The spin-Peierls phase of 
Fig. ^(e) orders into the exotic Neel phase of Fig. ||(c) 
where spins on tetrahedra A are collinear, while spins on 
tetrahedra B are coplanar. Note that, whereas the bond 
order we have considered is at q = 0, the concomitant 
spin order need not be. 

dN/da 



1/JS 



1/2JS 








JS 



2JS 



3JS 



4JS 



FIG. 5. Spectrum of spin waves of the elastic py- 
rochlore lattice for the collinear Neel state of Fig. |^(a) for 
J2 : Ji = 1 : 1 (solid line) and 3 : 4 (dashed line). In the 
undistorted magnet, J2 = Ji, half of the modes are at oj = 0. 
The dominant peak at ljq = 43 {Ji — J2) is due to a large 
number of string modes. Inset: Local mode on a hexagonal 
loop of parallel spins. For clarity, spins above and below the 
kagome plane (a [111] pyrochlore plane) are not shown. 

String modes. The extreme frustration of pyrochlore 
antiferromagnets is reflected in their spin- wave spectrum. 
A large degeneracy of the ground state leads to a macro- 
scopic number of zero-frequency modes. Although a 
Jahn-Teller distortion lifts the frustration and moves the 
zero modes to finite frequences, we find that the spin- 
wave spectrum contains a remnant of that degeneracy. 
Remarkably, a large fraction of spin waves can remain 
degenerate in the deformed magnet, appearing as a res- 
onance at the frequency luq — ASSJ. The effective dif- 
ference of exchange couplings dJ between satisfied and 
frustrated bonds Jij = J — 26*4(8^ • Sj) is due to the 
biquadratic exchange (^). The resonance is a local spin 
wave reminiscent of the weather-vane mode of the kagome 
antiferromagnet ]l5[ |. In the present context, the mode 
is confined to any line of parallel spins. Such lines are 
abundant in coUinearly ordered states: each spin has ex- 
actly two parallel neighbors. When successive spins on 
the string precess with a phase shift of tt, spins adja- 
cent to that line can remain unaffected because exchange 
field at their locations does not precess. (It can be shown 
that, despite the magnetoelastic coupling, harmonic spin 
waves and phonons are decoupled in any coUinearly or- 



dered magnet.) 

For various kinds of collinear ground states, the string 
mode can live on straight lines of parallel spins [Fig. ||(a)], 
spirals [Fig. ||(b)], irregular lines, or even closed loops. A 
mode living on a short loop is truly a local resonance: the 
local density of spin waves dN/ du contains a 5- function. 
Open strings are localized in only two directions and will 
therefore be seen in the local density of states as a strong 
van Hove singularity dA^/dcj c>c \uj — ujq\^^/^ (Fig. |^). The 
shape of the strings determines the form factor of the res- 
onance observable, e.g., by inelastic neutron scattering. 
A recently found resonance in ZnCr2Q4 has the form 



factor of a hexagonal string mode |16| (the shortest pos- 
sible loop. Fig. g, inset). 

Discussion. We have demonstrated that a large spin 
degeneracy and high symmetry of the pyrochlore network 
provide a robust way to relieve spin frustration: by means 
of a Jahn-Teller distortion of tetrahedra. A transition to 
a spin-Peierls state with a two-component order param- 
eter, which can be spatially nonuniform, may precede, 
or coincide with, a Neel transition. A unique feature of 
such a distorted antiferromagnet is a large number of de- 
generate spin waves confined to lines of parallel spins. 
Our theory has direct relevance to the observed mag- 
netic and structural phase transitions and local modes in 
pyrochlore antiferromagnets YMn2, ZnV2 04, MgV2 04, 
ZnCr204, and MgCr204. The extension of the above 
mechanism to other frustrated lattices such as the tri- 
angular and kagome is immediate for the case of Ising 
magnets; for Heisenberg symmetry the 5=1/2 kagome 
system with its high density of low energy states [^ ap- 
pears to be a promising candidate for future work. 
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For an introduction to f rustrated magnets, see R. Moess- 
ner, cond-mat/0010301; reviews of exact diagonalisations 
and experiments, respectively, are Ref. [^ and Q 
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